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Abstract 

We investigate the sensitivity of J /'^ production in proton-neutron and 
proton-deuteron collisions to charge symmetry violations in the parton dis- 
tributions as well as to nuclear corrections. It is found that the effects 
of charge symmetry violations in the quark distributions are confined to 
large xp and difficult to measure in experiments currently being planned. 
However, from proton-deuteron experiments it should be possible to isolate 
nuclear corrections to the gluon distribution. 
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1 Introduction 



The production of heavy quarkonia has been studied quite extensively in 
deep-inelastic scattering (see e.g. ref.|]l[|) and in hadron-hadron collisions 
(see e.g. ref. [||, |3|). In the latter, at sufficiently high energy, the reaction 
proceeds through the annihilation of either a quark or gluon in one hadron 
with an antiquark or gluon in the other hadron. In principle one therefore 
has a probe of the quark and gluon distributions. For the quarks this process 
merely supplements the information obtained in Drell-Yan and deep- inelastic 
scattering, but for the gluons this is one of our few direct probes. 

There have recently been suggestions that there may be a somewhat 
larger violation of charge symmetry in the minority valence quark distribu- 
tion in the nucleon than one might naively have expected. At intermediate 
X this effect could be as large as 5% |5|. It is important to try to test 
these ideas in Drell-Yan reactions Here we examine whether such an 
effect could induce a significant forward-backward asymmetry in the pro- 
duction of J/^ in proton- neutron collisions. As there has been no explicit 
calculation of a difference in the gluon distribution of a proton and neutron 
we also comment on this possibility. 

In proton-deuteron collisions the forward-backward asymmetry of J /ip 
production is sensitive not only to charge symmetry violations but also to 
nuclear modifications of the deuteron parton distributions. The latter could 
be investigated through experiments feasible at RHIC. 



2 J/ip production in proton- neutron collisions 



The production of charm-anticharm quark pairs proceeds in leading order in 
the strong coupling via gluon-gluon fusion or quark-antiquark annihilation. 
The corresponding cross sections for these subprocesses are Q: 
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where A = ^l- Am'i/AP. The mass of the charm quark and the produced 
cc pair is denoted by nic and M respectively. The strong coupling constant 
Os is calculated at M'^ using a QCD scale A = 0.177 GeF. Multiplying 
the cross sections of the QCD subprocesses with the parton distributions of 
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beam and target, fl' and /j, yields the cc production cross section 

In the center of mass frame Feynman Xi;' is defined as the fraction of beam 
momentum carried by the produced cc pair and s stands for the squared 
center of mass energy. The hght-cone momentum fractions of the active 
beam and target parton are 



1 4M^ 1 

The J/ij} production cross section is proportional to the charm produc- 
tion cross section, integrated over the invariant mass of the cc pair Q . The 
integration limits are the thresholds for cc pair production and DD produc- 
tion: 



dxpdM'^' 

The factor F specifies the fraction of events in which J /ip bound states are 
formed. Despite being a simple model, such a description is well known to 
describe many features of heavy quark production, including the dependence 
of the J production cross section on Xi;' and the beam energy |Q, |3[. Since 
it is of importance for our further discussion, we show in Fig. 1 the separate 
contributions to the J/ip production cross section from gluon fusion and 
quark-antiquark annihilation. We take rric = 1.5 GeV and use the parton 
distributions of ref.|^]. At small xp gluon fusion is by far the dominant 
mechanism, while for xp > 0.6 the qq annihilation takes over. 

To investigate charge symmetry violation we will now focus on J/ tp pro- 
duction in proton-neutron collisions. The corresponding production cross 
section involves the following combination of parton distributions 

^qq [Up{xi))Un (Xt) + Up{xb)Un{xt) + 

dp{xh)dn{xt) + dp{xb)dn{xt)] + 
^gg 9p{Xb)gn{xt) (7) 

Here are the quark distributions of the proton and neutron respectively 
and Qpjn the corresponding gluon distributions. 

In an isospin rotated world J /tp production in proton-neutron collisions 
becomes J /ip production in neutron-proton collisions, i.e. the role of beam 
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and target is interchanged. If charge symmetry were exact, the difference 
between the corresponding cross sections would vanish. Interchanging the 
role of beam and target is equivalent to a sign change in xp- Hence, the 
difference of the J /ip production cross sections at positive and negative xp 



dxp 



XF dxF 



(8) 



-Xp 



is driven by charge symmetry violations only. In detail, Ao"p„ contains the 
following combination of parton distribution functions: 

^qql { (Su{xh) - 6d{xb)) {d{xt) - u{xt)) + 

{5u{x}j) + 5d{xb)) {d{xt) + u{xt)) - 

{u{xb) - d{xb)) i6d_{xt) - 6u{xt)) - 

{u{xb) + d{xb)) {Sd{xt) + 6u{xt)) + 

[q < — > g]| + 

gg [S9{xb)9{xt) - 9{xb)Sg{xt)^ (9) 



cr. 



We expressed the neutron distributions through the proton ones, using the 
definitions 5d = dp — Un, 6u = Up — dn and Sg = gp — Qn, and dropping the 
index "p". 

First we will consider contributions to Aupn through charge symmetry 
violations in the valence distributions. The corresponding charge symmetry 
violating parts of the minority and majority distributions, 6d'" = d^ — 
and 5u^ = Up — d^ were extensively discussed in |@, ^ . We use the results of 
ref.Q which were obtained within the framework of the MIT bag model. The 
magnitude of 5d^ was found to be similar to that of du" . As d^ is generally 
much larger than the fractional change in d" is much greater. This can 
be easily understood, since one of the major sources of charge symmetry 
violation is the mass difference of the residual di-quark pair, when one quark 
of the nucleon is hit in a deep-inelastic scattering process. For the minority 
quark distribution the residual di-quark is uu in the proton and dd in the 
neutron. Therefore in the difference dp - Un the up-down mass difference 
enters twice. On the other hand, for the majority quark distribution, where 
the residual system is ud, both for the proton and the neutron, there is no 
contribution to charge symmetry breaking. 

To calculate A(jp„ we also need to know the flavor asymmetry of the 
quark sea, d — u, which enters in the first term of Eq.(|9|). We take a parame- 
terization from ref.Q: x{d{x) — u{x)) = Ax^'^{l — xY, where the normaliza- 
tion A is fixed through / dx{d{x) — u{x)) = 0.15. Such a parameterization is 
in good agreement with recently discovered violations of the Gottfried sum 
rule H. 
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We normalize Aapnixp) through the J/ip production cross section in 
proton-proton collisions and present in Fig. 2 results for the ratio Rpn^ = 
daiJ/^f" /dxp ^ center of mass energy ^/s = AOGeV. We find Rpi^ ~ 
—0.02 at large xp ~ 0.6 — 0.7. At smaller values of xp (say below 0.5) the 
cross section difference Acjp„ vanishes. This happens for two main reasons. 
At small xp gluon fusion yields the major contribution to J/ip production 
and quark-antiquark annihilation is of little relevance. Also 6d^ — du", which 
is much larger than Sd" + du" Q , enters Aupn in combination with the small 
flavor asymmetry, d — u. 

Charge symmetry breaking in the sea quark distributions has not been 
calculated yet. Nevertheless its influence on Aapn can be estimated. In the 
MIT bag model the sea quark distributions are dominated by contributions 
characterized through residual four-quark states, after one quark is hit in 
a deep-inelastic scattering process [10|. By analogy with the valence quark 
case let us assume that a major source of charge symmetry breaking is the 
mass difference of the residual spectator states. Scattering on an antiquark 
which carries the same flavor as the majority quarks leaves a uuud or dddu 
residual four-quark state in the case of a proton or neutron target respec- 
tively. In case of an antiquark of minority flavor a uudd residual state occurs 
in both cases. Since the up-down mass difference enters twice in the first 
case but is absent in the second, it seems reasonable to assume 6u ^ 6d. In 
the following we will neglect 6d. 

In Fig. 2 we also show Rpi^ for various values of 6u between O.Oln 
and O.lOn. Qualitatively we find minor changes to our former result where 
only charge symmetry breaking in the valence distributions was taken into 
account. We may therefore conclude that charge symmetry breaking in 
the quark distributions contributes to the difference of J/ip production in 
the forward and backward direction only at large values of xp (xp ^ 0.6). 
Unfortunately this kinematic region is not accessible at the moment. In 
current measurements, using the neutron beam facility at FNAL, one is 
restricted to the region xp > —0.1 |11]. 

However the news is not all bad. Since charge symmetry violations in 
quark distributions do not contribute to Rpd^ at small xp, it is an ideal 
place to look for charge symmetry violations in the gluon distributions. As 
we can see from Eq.(^, a (say) 1% charge symmetry violation in the gluon 

distribution, dg ~ 0.01(7, can lead to \Rpn^\ ^ 0.02. No predictions exist up 
to now for charge symmetry violations in gluon distributions. Since a major 
part of the glue in hadrons can be viewed as being radiatively generated 
from quarks, charge symmetry violations in the quark distributions may in 
principle induce similar effects in the gluon distributions. However, since 
for charge symmetry violation in the radiatively generated glue the small 
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combination Sd" + 6u'" is relevant, a large signal is not expected. If a large 
signal were seen it could only be attributed to charge symmetry violation in 
the non-perturbative glue, which would certainly be a surprise. 



3 J/i/j production in proton-deuteron collisions 

In high energy processes deuterons are often used as a convenient source of 
neutrons. Furthermore, in the near future proton-deutron collisions will be 
carried out at RHIC at center of mass energies between 50 and 375 GeV and 



xp > —0.5 [11]. Therefore, at first sight it seems to be a good idea to inves- 
tigate charge symmetry violations via J/'i/' production in proton-deuteron 
collisions. However, as we will demonstrate below, nuclear modifications of 
deuteron parton distributions, which are interesting in their own right, are 
most likely to overtake the effects resulting from charge symmetry violations. 

Let us assume that the parton distributions of the deuteron are related 
to those in the proton and neutron via 

qnix) = {I + eq{x)) {qp{x) + Qnix)), q = u,d, (10) 
9d{x) = {l + eg{x)){gp{x)+gn{x)). (11) 

In proton-deuteron collisions the difference Aapoixp) of the J/ip production 
cross sections measured in the forward and backward direction, involves the 
parton distribution functions: 

agg {1 + eg{xb)) [ {6u{xb) - 5d{xb)) {d{xt) - u{xt)) + 

{5u{xb) + 5d{xb)) {d{xt) + u{xt))] - 

\{l + eg{xt))[ {u{xb) - d{xb)) {5d{xt) - Su{xt)) + 

{u{xb) + d{xb)) iSd{xt) + 6u{xt))] + 

{eq{xt) - eq{xb)) iu{xb) + d{xb)) {u{xt) + d{xt)) + 

+ 



[q — ' ^]} 

I (1 + eg{xb)) 5g{xb) g{xt) - (1 + eg{xt)) g{xb) Sg{xt) + 
2 {(^g{xt) - eg{xb)) g(xb) g{xt) } 



(12) 



Clearly not all contributions to Aup^) arise from charge symmetry violation. 
There are also terms proportional to the difference of nuclear effects at Xb 
and xt- In the following we will estimate their size and show that they are 
most likely to dominate over contributions from charge symmetry violations. 

In Fig. 3 we show the light-cone momentum fractions Xb and xt at which 
the beam and target parton distributions are probed for different xp and 
for different center of mass energies ^/s. The invariant mass of the produced 
quark-antiquark pair is varied over the range 4m^ < < Am?j-) . We observe 
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that at large values of s the dependence of Xf/h on is rather small. While 
xt decreases to very small values with increasing xp, Xf, rises towards xp. 

Let us review the present knowledge of nuclear effects in deuteron dis- 
tribution functions at light-cone momentum fractions probed in J/-^ pro- 
duction. Nuclear effects in quark distributions have been studied a great 



deal in deep- inelastic scattering processes (see e.g. [|12|, |13[). At small val- 
ues of Bjorken x {x < 0.1) shadowing effects in the deuteron of about 
(1 — 4)% were predicted by many models and also recently observed (see 
[p^ , 15| and references therein). They suggest eq{x) ~ (—0.01) — (—0.04) 
for X ~ 0.05 — 0.001. At x ~ 0.1 shadowing disappears and the deep- 
inelastic scattering cross section for nuclear targets becomes slightly larger 
than the corresponding cross section for free nucleons. Such a behavior is 
also expected for deuterium, although small [|^. It suggests eq{x ~ 0.1) ^ 0. 
For X > 0.2 binding effects cause a decreases of eg < 0, while Fermi motion 
leads to a rise at x > 0.8 (see e.g. [17|). Nuclear effects in gluon distribu- 
tions are not so well established up to now. However momentum and baryon 
number conservation suggest that for x < 0.2 they exhibit a behavior similar 
to that of the quark distributions - i.e. eg ^ €q |16|. 

The preceeding discussion suggests that modifications of the deuteron 
parton distributions might easily be larger than the charge symmetry vio- 
lating effects. For example, at a center of mass energy = 80GeV and 
xf ~ 0.1 we have xt ~ 0.02 and Xb ~ 0.12. While xt is in the shadow- 
ing domain, xi, lies in the region where the deuteron parton distributions 
might be enhanced or are at least equal to the nucleon ones. The difference 
^^q/g ~ ^Q/g(^i) ~ ^q/gi^b) can therefore easily range from Ae^/^ —0.01 
to Aeq/g « -0.03 

In contributions to AapD which are proportional to charge symmetry 
violations one may neglect nuclear effects to a good approximation. Then 
nuclear effects enter AapD via the difference Ae^/g only. In Fig. 4 we show 
separately the contributions of gluon fusion and quark-antiquark annihila- 
tion to the ratio Rp^j^ = daijfii/)'^ /dxp ' different Ae^/^ , at a center of mass 
energy y/s = 80 GeV. For the charge symmetry violation in the valence dis- 
tributions we again use the results of ref.||5|, while the charge symmetry 
violation in the sea and gluon distributions are chosen to be zero. We find 
that at small values of xp nuclear modifications of the gluon distribution 
in deuterium dominate Rpj^ or equivalently Aapo- Therefore, if charge 
symmetry violation in the gluon distributions, which are accessible through 
Rpri^ , are small, nuclear modifications of the gluon distribution in deuterium 
can be investigated. Figure 4 demonstrates that, at large xp, nuclear effects 
may easily dominate charge symmetry violation in the valence distributions. 
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4 Conclusion 



We have discussed J/ijj production in proton-neutron and proton-deuteron 
collisions as a tool for investigating charge symmetry breaking and nuclear 
modifications of parton distributions. In proton-neutron collisions the dif- 
ference of the J/V' production cross section in the forward and backward 
direction, Acipn) is solely due to charge symmetry violations in the nucleon. 
Charge symmetry breaking in the valence and sea quark distributions affect 
the cross section difference significantly only at large xp 0.6. At small 
values oi xp ^ 0.1 a non vanishing result for A(7p„ would be entirely due 
to charge symmetry violations in the gluon distributions. Corresponding 
measurements should be possible using the neutron beam facility at FNAL. 

J /ip production in the forward and backward direction through proton- 
deuteron collisions will be possible at RHIC for a wide range of center of 
mass energies and Feynman xp- The corresponding cross section difference 
is, however, not only due to charge symmetry violations but also to nuclear 
modifications of the deuteron parton distributions. As a consequence nuclear 
modifications of the gluon distribution in the deuteron could be investigated 
as well. 

We would like to thank J. Moss for helpful discussions. 
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Figure 1: Ratio of the quark (solid) and the gluon (dashed) contribution 
to the J/ip production to their sum, for proton-proton colUsions at -y/s = 
AQGeV. 



Figure 2: Ratio Rpn at y/s = 40GeV. The sohd line includes charge 
symmetry violations in the valence quark distributions only. Different values 
for charge symmetry violations in the sea quark distributions are used: 6u = 
O.Olu (dashed), 5u = O.Obu (dotted) and 6u = O.lu (dot-dashed). 



Figure 3: Xi^-dependence of the light-cone momentum fractions (solid) 
and xt (dashed) for different center of mass energies. The shaded area 
indicates their dependence on the mass Am^ < < 4m|, of the final cc 
state. 



Figure 4: Gluon (dashed) and quark (dotted) contributions to the ratio 
^piD ~ SOGeV. For the nuclear modifications of the gluon and 

quark distributions in the deuteron Acg/g = —0.01, —0.02 are used. The 

solid line shows R^if in absence of nuclear effects. 
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